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Abstract 

Many  imaging  applications  require  the  acquisition  of  a 
time  series  of  images.  In  conventional  Fourier  transform- 
based  imaging  methods,  each  of  these  images  is  acquired 
independently.  Ai  a  result,  one  often  has  to  sacrifice  spatial 
resolution  for  temporal  resolution.  To  address  this  problem, 
this  paper  extends  a  generalized  series  imaging  method  so 
that  multiple  references  can  be  used  to  achieve  high  spa- 
tiotemporal  resolution.  Application  results  are  also  pre¬ 
sented  to  illustrate  its  effectiveness  for  high-resolution  dy¬ 
namic  imaging. 


1.  Introduction 


Magnetic  resonance  imaging  (MRI)  is  usually  imple¬ 
mented  as  a  Fourier  transform-based  technique.  As  such, 
the  imaging  equation  is  in  the  form  of  a  Fourier  integral. 
For  notational  convenience,  we  define  here  the  dynamic 
imaging  problem  as  the  acquisition  of  a  sequence  of  Q  im¬ 
ages,  denoted  as  A  (z),  A  (z),  •  •  • ,  Iq{x),  each  of  which  is  a 
snapshot  of  a  time-varying  image  function  I{x,t).  Conven¬ 
tionally,  Q  data  sets  in  A:-space,  di  (A:),  A  (A:),  ■  ■  ■  ,dQ{k), 
are  acquired  independently.  Assume  that  N  sample  data 
points  (or  encodings)  are  collected  for  each  of  these  data 
sets  at  A:  =  nAk,  —N 12  <  n  <  N 12.  The  spatial  res¬ 
olution  of  Iq{x)  is  limited  to  1/{N Ak),  whereas  the  tem¬ 
poral  resolution  is  NTn  with  Tn  being  the  repetition  time. 
Therefore,  high  spatial  resolution  requires  a  large  N ,  which 
means  poor  temporal  resolution. 

To  overcome  this  problem,  several  data-sharing  meth¬ 
ods  have  been  proposed  for  efficient  dynamic  imaging 
[1-6],  two  examples  of  which  are  the  Keyhole  and  RIGR 
(Reduced-encoding  Imaging  by  Generalized-series  Recon¬ 
struction)  techniques  [1-3].  This  paper  presents  an  exten¬ 
sion  on  the  RIGR  technique  so  that  multiple  references  can 
be  used  to  achieve  high  spatiotemporal  resolution. 


2.  Generalized  Series  Imaging 

2.1  The  generalized  series  model 

The  generalized  series  (GS)  model  is  a  general  mathe¬ 
matical  model  developed  for  constrained  image  reconstruc¬ 
tion.  In  this  model,  an  image  function  is  represented  as 

=  'Y^CnPn{x),  (1) 

n 

where  Pn{x)  are  basis  functions  adaptively  chosen  accord¬ 
ing  to  the  conditions  of  a  particular  application  problem. 
Of  special  interest  to  dynamic  imaging  is  the  class  of  basis 
functions  given  by 

^„(z)  =C(z)e*2-«^^L  (2) 

where  the  constraint  function  C(z)  is  normally  chosen  to 
absorb  available  a  priori  information.  Many  desirable  prop¬ 
erties  of  these  basis  functions  have  been  discussed  in  the 
context  of  constrained  spectroscopic  imaging  [7]  or  image 
reconstruction  in  general  [8].  The  optimality  of  these  basis 
functions  can  also  be  justified  from  the  principle  of  mini¬ 
mum  cross-entropy  [9].  After  the  basis  functions  Pn{x)  are 
chosen,  the  series  coefficients  can  be  determined  by  solving 
a  set  of  linear  equations.  More  specifically,  forcing  Igs{x) 
to  be  data-consistent  yeilds 

Ar/2-l 

dfrnAk)  =  Cndffm  —  n)  Ak]  (3) 

n=-N/2 

for  —N/2  <  m  <  N/2  —  1,  where  do(^)  is  the  Fourier 
transform  of  C(z).  Equation  (3)  can  be  solved  efficiently 
using,  for  example,  the  Levinson  algorithm  [10]. 

3.  The  RIGR  Technique 

The  basic  generalized  series  dynamic  imaging  scheme 
(or  RIGR)  is  characterized  by  the  acquisition  of  one  high- 
resolution  reference  image  and  a  series  of  low-resolution 
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dynamic  data  sets.  The  number  of  phase  encodings  for  the 
reference  image  is  chosen  to  satisfy  the  spatial  resolution 
requirement,  whereas  the  number  of  phase  encodings  for 
each  dynamic  data  set  is  chosen  to  give  the  desired  temporal 
resolution. 

Reconstruction  of  the  dynamic  images  is  accomplished 
by  directly  using  the  generalized  series  model.  Specifically, 
let  /ref(»)  be  the  high-resolution  reference  image.  The  basis 
functions  are  constructed  as 

=  |4ef(x)|e*“^^  (4) 

The  series  coefficients  c„  are  then  determined  as  described 
in  Section  2. 

A  notable  limitation  with  the  basic  RIGR  imaging 
method  is  due  to  the  fact  that  the  GS  model  has  as  many 
terms  as  the  number  of  dynamic  encodings  and,  as  a  result, 
dynamic  information  captured  by  the  GS  model  often  does 
not  have  the  same  spatial  resolution  as  the  reference  image. 
We  demonstrate  in  this  paper  that  this  problem  can  be  effec¬ 
tively  overcome  if  two  or  multiple  references  are  available. 

4.  RIGR  Imaging  with  Multiple  References 

In  some  dynamic  imaging  applications,  it  is  possible 
to  acquire  two  or  more  high-resolution  images  at  differ¬ 
ent  times  of  the  dynamic  imaging  process.  For  exam¬ 
ple,  in  contrast-enhanced  dynamic  imaging  of  breast  can¬ 
cer,  one  can  acquire  a  high-resolution  pre-contrast  image  as 
the  baseline  reference  and  another  high-resolution  reference 
image  when  the  contrast  agent  is  strongly  visible  in  the  im¬ 
age.  To  reconstruct  the  dynamic  images  between  the  two 
references,  the  TRIGR  (Two-reference  RIGR)  method  [11] 
uses  the  difference  image  between  the  two  references  to 
shape  the  basis  functions  of  the  GS  model.  More  specifi¬ 
cally,  let 

C{x)  =  |/2.ref(*)  -  -fl.ref(*)|  (5) 

and 

Adq{k)  =  dq{k)  -  dl,ref(fe)  (6) 

Dynamic  signal  changes  between  the  gth  data  set  and  the 
first  reference  image  are  expressed  in  terms  of  the  GS  model 
as 

M/2-1 

M,{x)=C{x)  Y.  (7) 

n=-M/2 

where  M  is  the  number  of  dynamic  encodings  for  each  dy¬ 
namic  data  set.  The  c„  can  then  be  determined  by  fitting 
the  model  to  the  difference  data  Adq{k)  given  in  Eq.  (6)  us¬ 
ing  the  algorithm  discussed  in  Section  2.  After  AIq{x)  is 
found,  Iq{x)  is  generated  by  simply  adding  AIq{x)  to  the 
reference  image.  That  is, 

Iq{x)  =  /l,ref(*)  +  AIq{x).  (8) 


This  processing  scheme  has  been  demonstrated  to  be  ca¬ 
pable  of  reconstructing  dynamic  signal  variations  at  high 
resolution  when  these  variations  are  of  an  additive  type  [11], 
such  as  the  case  with  dynamic  imaging  of  injected  contrast 
agents.  For  multiplicative  types  of  dynamic  signal  changes, 
as  is  the  case  with  diffusion-tensor  imaging,  this  scheme 
is  not  effective.  To  overcome  this  problem,  we  have  pro¬ 
posed  an  improved  TRIGR  reconstruction  method,  called 
TRIGR-h.  This  method  creates  two  sets  of  basis  functions 
based  on  /i,ref(®)  and  l2,mf{x),  respectively.  With  these  ba¬ 
sis  functions,  two  reconstructions  of  Iq{x)  are  obtained  as 

M/2-1 

Ii,q{x)  =  \h,M\  Y  (9) 

n=-M/2 

and 

M/2-1 

/2.,(x)  =  |/2.ref(x)|  Y  (10) 

n=-M/2 

The  series  coefficients  ci,„  and  C2,n  are  determined  by  fit¬ 
ting  Eqs.  (9)  and  (10)  to  dq{mAk)  =  dq{mAk)wm,  re¬ 
spectively,  where  Wm  is  chosen  according  to  Wm  =  0.54-1- 
0.46  cos(27rm/M),  for  —M/2  <  m  <  M/2.  This  filter 
function  serves  two  purposes:  (a)  to  reduce  ringing  artifacts 
due  to  model  truncation,  and  (b)  to  create  a  residual  signal 
Adq{mAk)  =  dq{mAk){l  —  Wm),  which  can  be  used  to 
pick  up  some  additional  dynamic  signal  changes  missed  by 
both/i,g(z)  and/2,g(z). 

Note  that  for  a  large  M,  both  Ii,q{x)  and  l2,q{x)  will 
reconstruct  Iq{x)  accurately.  This  is  not  the  case  when 
M  is  small.  In  fact,  in  this  case,  Ii,q{x)  will  be  different 
from  l2,q{x)  although  they  are  both  data-consistent.  One 
can  argue  that  l2,q{x)  —  Ii,q{x)  will  have  large  values  in 
areas  where  dynamic  signal  changes  are  not  well  repro¬ 
duced  by  Ii,q{x)  and/or  l2,q{x).  Therefore,  we  can  use 
\l2,q{X')  —  Ii^q{x)\  to  coustTuct  another  set  of  basis  func¬ 
tions  to  capture  these  dynamic  signal  changes.  Specifically, 
let 

M/2-1 

Alq{x)  =  \l2,q{x)-h,q{x)\  ^  (11) 

n=-M/2 

Fitting  this  model  to  the  residual  signal  Adq{mAk)  ob¬ 
tained  above  will  not  only  pick  up  additional  dynamic  signal 
changes  but  also  assure  the  final  reconstruction  to  be  con¬ 
sistent  with  dq{mAk)  for  —M/2  <  m  <  M /2. 

Specifically,  combining  Ii^q{x),  l2,q{x)  and  AIq{x)  as 
follows  will  yield  a  data-consistent  reconstruction  for  Iq{x)\ 

Iq{x)  =  aqll,q{x)  -f  (1  -  aq)l2,q{x)  +  AIq{x).  (12) 

The  weighting  coefficient  a,  can  be  determined  in  a  variety 
of  ways.  For  example,  assume  that  c/i,ref(a:)  and  c/2,ref(a:) 


are  acquired  at  i  =  ii^r  and  t2,i,  respectively,  and  that  dq{k) 
is  collected  at  t  =  tg  such  that  <  tg  <  i2,r-  We  may 
set  aq  =  (i2,r  —  tg)/{t2,i  —  ii,r)-  TWs  scheme  is  particu¬ 
larly  suitable  for  the  case  in  which  dynamic  signal  changes 
are  roughly  a  linear  function  of  time.  For  other  models  of 
dynamic  signal  changes,  the  ag  can  be  chosen  in  a  more  op¬ 
timal  fashion.  Nevertheless,  we  have  found  that  the  simple 
scheme  works  reasonably  well  for  a  number  of  problems, 
including  diffusion-tensor  imaging  and  cardiac  imaging. 

5.  Result 

Figure  1  illustrates  the  performance  of  the  proposed 
method.  In  this  experiment,  dynamic  signal  changes 
were  introduced  with  variable  diffusion  weightings  using 
a  diffusion-weighted  imaging  sequence.  The  images  in  the 
top  row  were  reconstructed  from  128  phase  encodings  (each 
with  128  points)  using  the  conventional  Fourier  method. 
The  diffusion  coefficient  map  was  obtained  by  fitting  a  sin¬ 
gle  exponential  curve  through  the  image  sequence  pixel  by 
pixel.  When  the  number  of  phase  encodings  is  reduced  to 
8,  the  Fourier  images  in  the  second  row  show  significant 
blurring  as  expected.  The  RIGR  method  using  the  same 
data  in  the  second  row  plus  a  high-resolution  reference  at 
the  beginning  of  the  sequence  produced  the  images  in  third 
row.  As  can  be  seen,  although  the  diffusion-weighted  im¬ 
ages  appear  to  be  high  resolution,  the  corresponding  dif¬ 
fusion  map  is  a  low-resolution  one,  which  means  that  the 
dynamic  signal  changes  were  not  produced  in  high  resolu¬ 
tion  in  the  diffusion-weighted  images.  Using  the  proposed 
method  and  the  same  data  in  the  second  row  plus  two  high- 
resolution  references,  one  at  the  beginning  and  the  other  at 
the  end  of  the  sequence,  we  obtained  the  images  in  the  last 
row.  These  images  demonstrate  that  the  method  is  able  to 
capture  dynamic  signal  changes  at  high  resolution  with  only 
a  small  number  of  encodings. 

6.  Conclusion 

An  improved  method  for  generalized  series  imaging  has 
been  discussed  in  this  paper.  This  method  enables  high- 
resolution  dynamic  images  to  be  reconstructed  from  a  re¬ 
duced  number  of  phase  encodings  plus  two  or  more  refer¬ 
ences.  We  have  applied  the  technique  to  several  applications 
which  demonstrate  that  a  4  to  8  fold  improvement  in  imag¬ 
ing  efficiency  over  the  traditional  Fourier  imaging  method 
can  be  obtained  without  a  significant  loss  of  image  qual¬ 
ity.  This  technique  should  prove  useful  for  improving  tem¬ 
poral  resolution  for  various  dynamic  imaging  experiments, 
such  as  dynamic  study  of  injected  contrast  agents,  diffusion- 
tensor  imaging,  and  cardiac  imaging. 
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Figure  1.  Diffusion-weighted  images  reconstructed  using:  (a)  the  Fourier  method  with  128  x  128 
encodings;  (b)  the  Fourier  method  with  8  x  128  encodings;  (c)  the  RiGR  method  using  the  same 
data  in  (b)  pius  a  high-resoiution  reference  at  the  beginning  of  the  sequence;  and  (d)  the  TRiGR+ 
method  using  the  same  data  in  (b)  pius  two  high-resoiution  references,  one  at  the  beginning  and  the 
other  at  the  end  of  the  sequence.  The  diffusion  maps  were  obtained  by  fitting  the  diffusion  curve 
pixei-by-pixei  through  the  corresponding  image  sequences  in  (a)-(d),  respectiveiy. 


